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Abstract 

A summary of how black holes grow in full, non-linear general relativity is presented. Specifically, 
a notion of dynamical horizons is introduced and expressions of fluxes of energy and angular 
momentum carried by gravitational waves across these horizons are obtained. Fluxes are local and 
the energy flux is positive. Change in the horizon area is related to these fluxes. The flux formulae 
also give rise to balance laws analogous to the ones obtained by Bondi and Sachs at null infinity 
and provide generalizations of the first and second laws of black hole mechanics. 

PACS numbers: 04.25.Dm, 04.70.Bw 



* Text of a plenary talk given at the conference on Mathematical Physics, General Relativity and Cosmology, 
held in honor of Professor Jerzy Plebanski in Mexico city, in September of 2002. To appear in the 
Proceedings. Results reported here were obtained jointly with Badri Krishnan and summarized very 
briefly in [1] . Details of proofs and extensions of these results will appear elsewhere. 
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I. INTRODUCTION 



Black holes are perhaps the most fascinating manifestations of the curvature of space and 
time predicted by general relativity. Properties of isolated black holes in equilibrium have 
been well-understood for quite some time. However, in Nature, black holes are rarely in 
equilibrium. They grow by swallowing stars and galactic debris as well as electromagnetic 
and gravitational radiation. For such dynamical black holes, the only known major result in 
exact general relativity has been a celebrated area theorem, proved by Stephen Hawking in 
the early seventies: if matter satisfies the dominant energy condition, the area of the black 
hole event horizon can never decrease. This theorem has been extremely influential because 
of its similarity with the second law of thermodynamics. However, it is a 'qualitative' result; 
it does not provide an explicit formula for the amount by which the area increases in any 
given physical situation. One might hope that the change in area is related, in a direct 
manner, to the flux of matter fields and gravitational radiation falling into the black hole. Is 
this in fact the case? If so, the formula describing this dynamical evolution of the black hole 
would give us a 'finite' generalization of the first law of black hole mechanics: The standard 
first law, 8E = (n/^G)5a + 0,5 J, relates the infinitesimal change 5a in the black hole area 
due to the infinitesimal influx 5E of energy and angular momentum 5 J as the black hole 
makes a transition from one equilibrium state to a nearby one, while the exact evolution 
law would provide its 'integral version', relating equilibrium configurations which are far 
removed from one another. 

From a general, physical viewpoint, these expectations seem quite reasonable. Why, then, 
had this question remained unresolved for three decades? The reason is that when one starts 
thinking of possible strategies to carry out these generalizations, one immediately encounters 
severe difficulties. To begin with, to carry out this program one would need a precise notion 
of the gravitational energy flux falling into the black hole. Now, as is well known, in full 
general relativity, there is no gauge invariant, quasi- local notion of gravitation radiation. The 
standard notion refers to null infinity, where one can exploit the weakness of curvature to 
introduce the notion of asymptotic translations, define Bondi 4-momentum as the generator 
of these translations, calculate energy fluxes, and prove the 'balance law' relating the change 
in the Bondi 4-momentum with the momentum flux across portions of null infinity. Even 
this structure at null infinity is highly non-trivial. Indeed there was considerable confusion 
about physical reality of gravitational waves well into the late fifties because it was difficult 
to disentangle coordinate effects from physical ones. Even Arthur Eddington who, according 
to a popular legend, was 'one of the only three wise men' to understand general relativity 
from the very early days, is said not to have believed in the reality of gravitational waves. 
Apparently, he referred to them as pure coordinate effects which 'travelled at the speed of 
thought'! Therefore, the introduction of the Bondi framework in the sixties was hailed as 
a major breakthrough: it established, once and for all, the physical reality of gravitational 
waves. Bondi famously said: They are real; they carry energy; you can heat water with 
them! 

However, to address the issues discussed above pertaining to black holes, one does not 
have the luxury of working in the asymptotic region; one must face the highly curved ge- 
ometry near black holes. Since this is the strong field, highly non-linear regime of general 
relativity, it seems hopeless to single out a translation group unambiguously. What would 
the notion of energy even mean without the usual simplifications at infinity? There do 
exist formulas for the amount of gravitational energy contained in a given region. But typi- 
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cally they use pseudo-tensors and are coordinate and/or frame dependent. Consequently, in 
strong curvature regions, these expressions fail to be gauge invariant whence, from a physi- 
cal perspective, they are simply not meaningful. Thus, even a broad conceptual framework 
or paradigm was not available within which one could hope to establish a formula relating 
the change in area to the flux of energy and angular momentum falling into the black hole. 
The issues had remained unresolved because it appeared that one has to simultaneously de- 
velop the conceptual framework which is to provide a natural home for the required notions 
in the strong field regime of a black hole, and manipulate field equations with their full 
non-linearity to get explicit expressions for energy and angular momentum fluxes. 

At first these challenges seem formidable. However, it turns out that an appropriate 
paradigm is in fact suggested by the strategy used in numerical simulations of black hole 
formation and merger. Using it, a program was initiated in collaboration with Badri Kr- 
ishnan [1, 2] a few months before the Plebanskifest and has been further developed by the 
two of us as well as by Ivan Booth and Stephen Fairhurst since then. The final results 
are surprising because one can introduce the necessary notions of energy and angular mo- 
mentum fluxes in the strong field, fully non-linear regime. But the methods used are all 
well-established and quite conservative; there is nothing here that was not available in the 
seventies. 

The purpose of this contribution is to provide a bird's eye view of the present situation 
and of prospects for near future. A detailed and more comprehensive treatment will appear 
elsewhere. 

II. CONCEPTUAL PARADIGM 
A. The idea 

Our first task is to give a precise definition of the black hole surface whose area is to 
increase during evolution. 

Heuristically, one thinks of black holes as regions of space-time from which no signal can 
escape to infinity. In mathematical general relativity, this idea is captured in the notion of 
event horizons. Standard space-times — such as the Kerr family — that we physically think 
of as containing black holes, all have event horizons. However, the notion of event horizons 
is extremely non-local and teleological: it is the future boundary of the past of future null 
infinity. Consequently, one can use event horizons to detect the presence of a black hole 
only if one has access to the full space-time metric, all the way to the infinite future. This 
extreme non-locality is dramatically illustrated by the fact that there may well be an event 
horizon developing in your room as you read this page because a million years from now, 
there may be a gravitational collapse in a nearby region in our galaxy! 

Because of these features, generally, the notion of an event horizon is not very useful in 
practice. A striking example is provided by numerical simulations of a gravitational collapse 
leading to the formation of a black hole (or, of binary black holes which merge). Here, 
one is interested in evolving suitable initial data sets and one needs to know, at each time 
step, if a black hole has formed and, if so, where it is. The teleological nature of the event 
horizon makes it totally unsuitable to detect black holes during these evolutions. One needs 
a notion that can sense the presence of the black hole and narrow down its approximate 
location quasi-locally, e.g., using only the initial data that can be accessed at each instant of 
time. Marginally trapped surfaces and apparent horizons provide such notions. A marginally 
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trapped surface S is a 2-dimensional sub-manifold of space-time A4, topologically S* 2 , such 
that the expansion of one null normal to it, say £ a , is everywhere zero and that of the other 
null normal, n a , is negative. The second condition merely says that n a is the inward pointing 
null normal. The non-trivial feature of S is that the expansion of the other null normal, £ a , 
is zero rather than positive. Thus, none of the light rays emerging from any point on S are 
directed towards the 'asymptotic' or 'outside' region. Apparent horizons are associated with 
(partial) Cauchy surfaces M. Given M, the apparent horizon S is the outermost marginally 
trapped surface lying in M. In numerical simulations then, one keeps track of black holes by 
monitoring the behavior of apparent horizons that emerge during the evolution. Assuming 
cosmic censorship, once an apparent horizon develops, there is an event horizon which lies 
outside it. 

Thus, in numerical relativity, one typically has a foliation of the given region of the 
space-time M. by partial Cauchy surfaces M, each equipped with an apparent horizon. The 
apparent horizons can 'jump' discontinuously. For example, in the black hole coalescence 
problem, there are two distinct apparent horizons up to a certain time step and then there is 
a sudden jump to a single connected apparent horizon. However, in practice, these discrete 
jumps happen only at a few places during the course of numerical simulation. Here, we will 
be concerned with those evolutionary epochs during which the 'stack of apparent horizons' 
formed by evolution span a continuous world tube. Along these apparent horizon world- 
tubes, the area increases and our aim is to present formulas relating this increase with the 
influx of energy and angular momentum. 

Before proceeding to give precise definitions, however, let me emphasize that the above 
considerations suggested by numerical relativity serve only as motivation. In particular, we 
will not need a foliation of space-time M. by partial Cauchy surfaces M. The object of 
direct physical interest — a dynamical horizon — can be located using only the space-time 
geometry, although in applications to numerical relativity, it will typically arise as the world 
tube of apparent horizons. 

B. Definition and methodology 

Let us then begin with the precise definition. 
Definition: A smooth, three-dimensional, space-like sub-manifold H in a space-time M. is 
said to be a dynamical horizon if it can be foliated by a family of two-spheres such that, 
on each leaf S, the expansion 9^ of a null normal £ a vanishes and the expansion 0( n ) of the 
other null normal n a is strictly negative. 

Thus, a dynamical horizon if is a 3-manifold which is foliated by marginally trapped 
2-spheres. Note first that, in contrast to event horizons, dynamical horizons can be located 
quasi-locally; knowledge of the full space-time is not required. Thus, for example, you can 
rest assured that no dynamical horizon has developed in the room you are sitting in ever 
since it was built! On the other hand, while in asymptotically flat space-times black holes are 
characterized by event horizons, there is no one-to-one correspondence between black holes 
and dynamical horizons. It follows from properties of trapped surfaces that, assuming cosmic 
censorship, dynamical horizons must lie inside an event horizon. However, in the interior 
of an expanding event horizon, there may be many dynamical horizons. 1 Nonetheless, the 



At this still preliminary stage, there is essentially no control on how many dynamical horizons there 
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framework is likely to have powerful applications to black hole physics because its results 
apply to each and every one of these dynamical horizons. 

Apart from the requirement that H be foliated by marginally trapped surfaces, the def- 
inition contains only two conditions. The first asks that H be space-like. This property 
is implied by a stronger but physically reasonable restriction that the 'inward' derivative 
C n 9m of 0m be negative and the flux of energy across H be non-zero. The second condition 
is that the leaves be topologically S 2 . This can be replaced by the weaker condition that 
they be compact. One can then show that the topology of S is necessarily S* 2 if the flux of 
matter or gravitational energy across H is non-zero. 2 Thus, the conditions imposed in the 
definition appear to be minimal. If these fluxes were to vanish identically, H would become 
isolated and replaced by a null, non-expanding horizon [5]. 

Dynamical horizons are closely related to Hayward's [3] future, outward trapping horizons 
H' . These H' are 3- manifolds, foliated by compact 2- manifolds S with 9(e) = 0, 9( n ) < 0, 
and C n 9(e) < 0. Assuming that the null energy condition holds, one can show that H is 
either space-like or null. If it is space-like, it is a dynamical horizon, while if it is null it is 
a non-expanding horizon, studied extensively in [5]. This notion is especially useful in the 
study of 'weakly-dynamical' horizons [4] which can be viewed as perturbations of isolated 
horizons [5-8]. The main difference from the notion of dynamical horizons used in this 
report is that while the definition of trapping horizons imposes a condition on the derivative 
of 9(e) off H, dynamical horizons refer only to geometric quantities which are intrinsically 
defined on H . But this difference is not physically significant because, in cases of interest, 
the additional condition would probably be satisfied and dynamical horizons will be future, 
outer trapping horizons. Rather, the significant difference with respect to Hayward's work 
lies in the way we analyze consequences of these conditions and in the results we obtain. 
While Hayward's framework is based on a 2+2 decomposition, ours will be based on the 
ADM 3+1 decomposition. We use different parts of Einstein's equations, our discussion 
includes angular momentum, our flux formulae are new and our generalization of black hole 
mechanics is different. 

Let us begin by fixing notation. Let r a be the unit time-like normal to H and denote by 
V the space-time derivative operator. The metric and extrinsic curvature of H are denoted 
by q a b and K ab := q a c Q.b d V c T d respectively; D is the derivative operator on H compatible with 
q a b and TZ^ its Ricci tensor. If H admits more than one foliation by marginally trapped 
surfaces, we will simply fix any one of these and work with it. Thus, our results will apply to 
all such foliations. Leaves of this foliation are called cross-sections of H. The unit space-like 



can be inside an evolving event horizon, but the expectation is that there will be several distinct ones. 
Stationary event horizons, by contrast, do not admit any dynamical horizons because, for simplicity, we 
have tailored our definition to cases where the area increases monotonically. In the closely related notion of 
outer trapping horizons, introduced by Hayward and discussed below, the area either increases or remains 
constant. So, the event horizon of a stationary space-time is an example of a trapping horizon. 
2 This may seem surprising at first. As Carlo Rovelli asked after my talk, in the problem of black hole 
merger, initially one has two distinct horizons and finally there is only one. So, the topology changes. 
How can one reconcile this with the universality of topology of HI Note first that to address this issue 
we should examine apparent — not event — horizons. Secondly, and more importantly, by definition H 
can be realized only by continuous segments of the world tubes of apparent horizons. Before the merger, 
each continuous segment does have topology S 2 x R and after the merger the final world tube has the 
same topology. In between the apparent horizon simply jumps and there is no dynamical horizon. 
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vector orthogonal to S and tangent to H is denoted by r a . Quantities intrinsic to S will be 
generally written with a tilde. Thus, the two- metric on S is q ab , the extrinsic curvature of 
S C H is K ab := q a c q b d D c 1 r ( i, the derivative operator on (S, q ab ) is D and its Ricci tensor is 
lZ a b- Finally, we will fix the rescaling freedom in the choice of null normals via £ a := r a + r a 
and n a := f a — r a . 

We first note an immediate consequence of the definition. Since 9^ = and 6*( n ) < 0, it 
follows that 

K = q ab D a r b = l - q ab V a {£ b - n b ) > 0. 

Hence the area a$ of S increases monotonically along r a . Thus the second law of black hole 
mechanics holds on H. Our first task is to obtain an explicit expression for the change of 
area. 

Our main analysis is based on the fact that, since if is a space-like surface, the Cauchy 
data (q ab , K ab ) on H must satisfy the usual scalar and vector constraints 

H s ■= n + K 2 - K ah K ab = lQirGT ab ? a ? b (2.1) 
H$ := D b (K ab - Kq ah ) = 8nGT bc r c q\ . (2.2) 

We will often fix two cross-sections Si and 5*2 and focus our attention on a portion AH C H 
which is bounded by them. 



III. ENERGY FLUXES AND AREA BALANCE 

Let us now turn to the task of relating the change in area to the flux of energy across H . 

As is usual in general relativity, the notion of energy is tied to a choice of a vector field. 
The definition of a dynamical horizon provides a preferred direction field; that along i a . To 
extract a vector field, we need to fix the proportionality factor, or the 'lapse' N, let us first 
introduce the area radius R, a function which is constant on each S and satisfies as = 4:7iR 2 . 
Since we already know that area is monotonically increasing, R is a good coordinate on H. 
Now, the 3-volume d 3 V on H can be decomposed as d 3 V = \8R\~ l d 2 VdR. Therefore, as we 
will see, our calculations will simplify if we choose N R = \dR\. Let us begin with this simple 
choice, obtain an expression for the change in area and then generalize the result to include 
a more general family of lapses. 

Fix two cross sections Si and S2 of H and denote by AH the portion of H they bound. 
We are interested in calculating the flux of energy associated with = Nf>£ a across AH. 

Denote the flux of matter energy across AH by Tm^'- 

?L R) := / T ab ? a e {R) d 3 V. (3.1) 

J AH 

By taking the appropriate combination of (2.1) and (2.2) we obtain 

?L R) = J^J ah N k i H s + 2W} <?V . (3.2) 

Since H is foliated by two-spheres, we can perform a 2 + 1 split of the various quantities on 
H. Using the Gauss Codazzi relation we rewrite 1Z in terms of quantities on S: 

TZ = TZ + K 2 ~ K ab K ab + 2D a a a (3.3) 
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where a a = r b D b r a — r a D b r b . Next, the fact that the expansion 6^ of £ a vanishes leads to 
the relation 

K + K = K ab r a r b . (3.4) 
Using (3.3) and (3.4) in (3.2) and simplifying, we obtain the result 

/ N R Ud 3 V = lQirG [ T ab r a ^ R) d 3 V+ [ N R { \a\ 2 + 2\(\ 2 } d 3 V (3.5) 
Jah Jah Jah 

where |<r| 2 = o ab o ab with a ab being the shear of £ a , and |(| 2 = ( a ( a with ( a := q ab r c V c l b) 
both a ab and ( a are tensors intrinsic to S. To simplify the left side of this equation, recall 
that the volume element d 3 V on H can be written as d 3 V = N R l dRd 2 V where d 2 V is the 
area element on S. Using the Gauss-Bonnet theorem, the integral of N R TZ can then be 
written as r 

[ N R TZd 3 V = [ 2 dR[ (f nd 2 v) = 8n{R 2 -Ri). (3.6) 
Jah Jr.! \Js J 

(It is this manipulation that dictated our choice of N R .) Substituting this result in (3.5) we 
finally obtain 

(§-§)= L T ^> d ' v + L Nr {h ' 2 + 2|c|2} sv - (3J) 

This is the first key result we were looking for. Let us now interpret the various terms 
appearing in this equation. The left side gives us the change in the horizon 'radius' caused 
by the dynamical process under consideration. The first integral on the right side of this 
equation is the flux T m of matter energy associated with the vector field £" R y Since £^ 
is null and r time-like, if T ab satisfies, say, the dominant energy condition, this quantity is 
guaranteed to be non-negative. Since the second term is purely geometrical and emerged 
as the 'companion' of the matter term, it is tempting to interpret it as the flux of 
£^-energy in the gravitational radiation: 3 

^=^/^MM 2 + 2|CI^V. (3.8) 

Is this proposal physically viable? We will now argue that the answer is in the affirmative 
in the sense that it passes all the 'text-book' tests one uses to demonstrate the viability of 
the Bondi flux formula at null infinity. 

First, since we did not have to introduce any structure, such as coordinates or tetrads, 
which is auxiliary to the problem, the expression is obviously 'gauge invariant'. Second, 
the energy flux is manifestly non-negative. Third, all fields used in it are local; we did not 
have to perform, e.g., a radial integration to define any of them. Fourth, the expression 
vanishes in the spherically symmetric case: Since the only spherically symmetric vector field 
and trace-free, second rank tensor field on a 2-sphere are the zero fields, if the Cauchy data 



3 While the presence of the shear term \a\ 2 in the flux formula (3.8) is natural from one's expectations 
based on the weak field limit, the term \(\ 2 is surprising. Booth and Fairhurst have shown that this term 
is of two orders higher than the shear term in the weak field expansion. Thus, it captures some genuinely 
non-linear, strong field physics which is yet to be understood fully. 
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(q a b, K a t,) and the foliation on H are spherically symmetric, a a b = and ( a = 0. Next, 
one might be concerned that the flux may not vanish in stationary space-times. Even in 
the Schwarzschild space-time, could one not construct a clever, non-spherical dynamical 
horizon HI If one could, the area law (3.8) would hold and then we would be led to an 
absurd conclusion that there is flux of gravitational energy across this H\ Even if this is not 
possible in the Schwarzschild space-time, could it not happen in a more general stationary 
space-time? If it can, the 'gravitational energy-flux' interpretation would not be viable. 
Now, since H is foliated by marginally trapped surfaces, it follows from general results that 
it must lie inside the event horizon. Using the fact that there is a Killing field in that 
region, one can show that there are no dynamical horizons in this interior region, whence 
the concern is unfounded. Thus, the expression on the right side of (3.8) shares with the 
Bondi-Sachs energy flux at null infinity all its key properties. We will therefore interpret 
it as the £(^)-energy flux of carried by gravitational waves. Recently, Booth and Fairhurst 
[4] have verified that on 'weakly-dynamical' horizons, the expression reduces to the familiar 
one from perturbation theory. They have also shown that this formula can be derived 
from a Hamiltonian framework where H is treated as the inner boundary of the space-time 
region of interest. These results provide considerable further support for our interpretation. 
Nonetheless, it is important to continue to think of new criteria and make sure that (3.8) 
passes these tests. 

Remark: The emergence of a precise formula for the flux of energy across AH is very 
surprising. What would happen if we repeat the above procedure for a general space-like 
surface HI The analog of the flux term would be much more complicated and fail to be 
positive definite. This happens even if we assume that H is foliated by strictly — rather than 
marginally — trapped surfaces S, i.e. if we replaced the condition am = by am < 0. Thus, 
there is no satisfactory candidate for the flux formula across H. To summarize, although 
the calculation is straightforward, it crucially depends on subtle cancellations which occur 
precisely because if is a dynamical horizon. 

To conclude this section, let us discuss the possibility of choosing more general lapse 
functions. In the above calculation, we needed a specific form of the lapse to cast Nro1 3 V 
into the form d 2 VdR. This suggests that we use more general functions r which are constant 
on each leaf S of the foliation and set N r = \dr\. If we use a different radial coordinate r', 
then the lapse is rescaled according to the relation 

N r i = —N r . (3.9) 
dr 

Thus, although the lapse itself will in general be a function of all three coordinates on H, 
the relative factor between any two permissible lapses can be a function only of r. Recall 
that, on an isolated horizon, physical fields are time independent and null normals — which 
play the role of N r £ a there — can be rescaled by a positive constant [5, 7]. In the present 
case, the horizon fields are 'dynamical', i.e., r-dependent, and the rescaling freedom is by 
a positive function of r. Thus, the freedom in the choice of lapse is just what one would 
expect. 

Given a lapse N r , following the terminology used in the isolated horizon framework, 
the resulting vector fields £^ := N r £ a will be called permissible. By repeating the above 
calculation, it is easy to arrive at a generalization of (3.8) for any permissible vector field: 

= L T " 9 '& dH ' + L n * { w? + 2|c|2} d * v • {3M) 
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where the constants r\ and r2 are values the function r assumes on the fixed cross-sections 
Si and This generalization of (3.7) will be useful in section V. 



IV. ANGULAR MOMENTUM 



To obtain the integral version of the full first law, we need the notion of angular mo- 
mentum and angular momentum flux. It turns out that the angular momentum analysis 
is rather straightforward and is, in fact, applicable to an arbitrary space-like hypersurface. 
Fix any vector field ip a on H which is tangential to the cross-sections of H . Contract tp a 
with both sides of (2.2). Integrate the resulting equation over the region AH C H, perform 
an integration by parts and use the identity £ ^qw, = 2D( a ip b ) to obtain 

= L( T ° brV + ^G p ° bc ^) d3y (44) 

where P ab := K ab — Kq ab . It is natural to identify the surface integrals with the generalized 
angular momentum associated with those surfaces and set: 

J s ] = f s K ^ a ? b d 2 V (4.2) 

where we have chosen the overall sign to ensure compatibility with conventions normally 
used in the asymptotically flat context. The term 'generalized' emphasizes the fact that the 
vector field tp a need not be an axial Killing field even on S; it only has to be tangential to 
our cross-sections. 

The fluxes of this angular momentum due to matter fields and gravitational waves are, 
respectively, 



Jg> = - / T ab ? a <p b d 3 V, (4.3) 

J AH 

and we get the balance equation 

_ = jM + jW) . (4.5) 

As expected, if tp a is a Killing vector of the three-metric q a b, then the gravitational flux 
vanishes: J^f 1 = 0. For the discussion of the integral version of the first law, it is convenient 
to introduce the angular momentum current 

f := -K abV a r b (4.6) 

so that the angular momentum formula becomes 



J. 



if) 



s 
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V. FINITE VERSION OF THE FIRST LAW 



Let us now combine the results of sections III and IV to obtain the physical process 
version of the first law for H and a mass formula for an arbitrary cross-section of H . 

To begin with, let us ignore angular momentum and consider the vector field £(m of 
section III. Denote by E^ R ~> the £(R)-energy of cross-sections S of H . While we do not yet 
have the explicit expression for it, it is natural to assume that, because of the influx of matter 
and gravitational energy, E^ R ) will change by an amount A = jffi + F ( g R) as we move 
from one cross-section to another. Then, the infinitesimal form of (3.7), dR/2G = dE^ R \ 
suggests that we define effective surface gravity associated with £? R -. as := 1/2R so 

that the infinitesimal expression is recast into the familiar form (kn/87iG)da = dE^ R ) where 
a is the area of a generic cross-section. For a general choice of the radial coordinate r, (3.10) 
yields a generalized first law: 

da = dE^) (5.1) 

87rG 

provided we define the effective surface gravity k r of £^ by 

dr dr 
K= m k R where = (5.2) 

Note that this rescaling freedom in surface gravity is completely analogous to the rescaling 
freedom which exists for Killing horizons, or, more generally, isolated horizons [5, 7]. The 
new feature in the present case is that we have the freedom to rescale the £ a and the surface 
gravity by a function of the radius R rather than by a constant. This is just what one 
would expect in a dynamical situation since R plays the role of 'time' along H. Finally, note 
that the differentials appearing in (5.1) are actual, physical variations along the dynamical 
horizon due to an infinitesimal change in r and are not variations in phase space as in the 
formulations [5, 7, 10] of the first law based on Killing or isolated horizons. Thus, (5.1) is 
a physical version of the first law, whence (3.10) is the finite version of the first law in the 
absence of rotation. 

Next, let us include rotation. Pick a vector field (f a on H such that (f a is tangent to the 
cross-sections of H, has closed orbits and has afline length 2tt. (At this point, ip a need not 
be a Killing vector of q a b.) Consider time evolution vector fields t a which are of the form 
t a = N r £ a — Vtip a where N r is a permissible lapse associated with a radial function r and 
f2 an arbitrary function of r. 4 (On an isolated horizon, the analogs of these two fields are 
constant.) Evaluate the quantity J AH T ab r a t b d 3 V using (4.1) and (3.7): 




- I vtfd 2 v- [ 2 dn I j v d 2 v\ = 

[ T ab r a t b d 3 V + -J— [ A r (H 2 + 2|C| 2 ) d 3 V--^— [ QP ab C ^d'V (5.3) 

J AH 1D7TCt J ah 107TG J ah 
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Constancy of £1 on cross-sections implies rigid rotation, although the frequency of rotation is allowed 
to change in 'time'. After my talk, Alberto Garcia raised the interesting issue of allowing differential 
rotations. This can be done by letting Q be a general function on H . In this case, one can still obtain 
the integral version of the first law but, as one would expect, if fl has angular dependence, one can not 
recover the familiar infinitesimal form of the first law. 
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This is our finite version of the familiar first laws of the isolated horizon framework [5, 7]. 
For, if we now restrict ourselves to infinitesimal AH, the three terms in the curly brackets 
combine to give d(QJ) — JdVL and we obtain 

(IT' n' 

— + ttdJ= ~-^da + ttdJ = dE\ (5.4) 

This equation is just the familiar first law but now in the setting of dynamical horizons. 
Since the differentials in this equation are variations along H, this can be viewed as a 
physical process version of the first law. Note that for each allowed choice of lapse N r , 
angular velocity f2(r) and vector field ip a on H, we obtain a 'permissible' time evolution 
vector field t a = N r £ a — Vt(p a and a corresponding first law. This situation is very similar 
to what happens in the isolated horizon framework [5, 7] where we obtain a first law for 
each permissible time translation on the horizon. Again, the generalization from that time 
independent situation consists of allowing the lapse and the angular velocity to become 
r-dependent, i.e., 'dynamical'. 

For every allowed choice of (N r , £l(r), y? a ), we can integrate (5.4) on H to obtain a formula 
for E l on any cross section but, in general, the result may not be expressible just in terms 
of geometric quantities defined locally on that cross-section. However, in some physically 
interesting cases, the expression is local. For example, in the case of spherical symmetry, it is 
natural to choose Q = and R as the radial coordinate in which case we obtain E t = R/2G. 
This is just the irreducible (or Hawking) mass of the cross-section. Even in this simple case, 
(5.3) provides a useful balance law, with clear-cut interpretation. Physically, perhaps the 
most interesting case is the one in which q ab is only axi-symmetric with ip a as its axial Killing 
vector. In this case we can naturally apply, at each cross-section S of H, the strategy used 
in the isolated horizon framework to select a preferred t a : Calculate the angular momentum 
J defined by the axial Killing field ip, choose the radial coordinate r (or equivalently, the 
lapse N r ) such that 

k r = UR) := , = (5.5) 

and choose f2 such that 

9C T 

n = n o (R):= (5.6) 

This functional dependence of k r on R and J is exactly that of the Kerr family. With this 
choice of N r and Q, the energy Eg is given by the well known Smarr formula 



E t °-2( k ° a | n J) + A ° 2J2 (57) 

Thus, as a function of its angular momentum and area, each cross-section is assigned simply 
that mass which it would have in the Kerr family. This may seem like a reasonable but rather 
trivial strategy. The non-triviality lies in two facts. First, with this choice, there is still a 
balance equation in which the flux of gravitational energy is local and positive definite 
(see (5.3)). (The gravitational angular momentum flux which, in general, has indeterminate 
sign vanishes due to axi-symmetry.) Second, as mentioned in section III, Booth and Fairhurst 
have recently shown that this expression of the dynamical horizon energy emerges from a 
systematic Hamiltonian framework on M. where H is treated as an inner boundary. 
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Motivated by the isolated horizon framework, we will refer to this canonical E to as the 
mass associated with cross-sections S of H and denote it simply by M. Thus, among the 
infinitely many first laws (5.4), there is a canonical one: 



Note that the mass and angular momentum depend only on geometrical fields on each cross 
section and, furthermore, the dependence is local. Yet, thanks to the constraint part of 
Einstein's equations, changes in mass over finite regions AH of H can be related to the 
expected matter fluxes and to the flux of gravitational radiation which is local and positive. 

I will conclude this section with a conceptual subtlety, emphasized by Stephen Fairhurst 
in recent conferences. The first law (5.3) discussed here is really a conservation law and 
as such it is a finite or integral version of the first law of black hole mechanics. In the 
discussion of the first law in its standard, infinitesimal form, one explicitly or implicitly 
considers transitions between an equilibrium state and a nearby equilibrium state. Concep- 
tually, this is the same setting as in laws of equilibrium thermodynamics. In particular, in 
infinitesimal processes involving black holes, the change in surface gravity can be ignored 
just as the change in the temperature is ignored in the first law dE = TdS + W of ther- 
modynamics. During fully non-equilibrium thermodynamical processes, by contrast, the 
system does not have time to come to equilibrium and there is no canonical notion of its 
temperature. Similarly, in the case of dynamical horizons, we only have a notion of 'effec- 
tive' or 'average' surface gravity; in striking contrast to what happens on isolated horizons 
which describe equilibrium configurations, K r does not have the geometrical interpretation 
of surface gravity. However, if one considers 'weakly-dynamical' horizons and regards them 
as perturbations of isolated horizons, there is a geometrical notion of surface gravity (of 
which K r is a 2-sphere average). In this situation, the geometrical surface gravity appears 
to be a good analog of the temperature, the idea being that system is evolving slowly so 
that it can reach approximate equilibrium in spite of time dependence. In these situations, 
the dynamical first law (5.3) can be simplified by keeping terms only up to second order in 
perturbation [4] and that expression can be regarded as the integral version of the first law 
of black hole thermodynamics. 

VI. DISCUSSION 

I will conclude the report by pointing out some important open problems whose solu- 
tions will add very significantly to our knowledge of dynamical horizons and suggest some 
applications of this framework. 

i) An important open problem is to obtain a complete characterization of the 'initial data' 
on dynamical horizons. Can one characterize the solutions to the constraint equations such 
that (H, q a b, K a b) is a dynamical horizon? Note that this would, in particular, provide a 
complete control on the geometry of the world tube of apparent horizons that will emerge in 
all possible numerical simulations! One can further ask: Can one isolate the freely specifiable 
data in a useful way? Are these naturally related to the freely specifiable data on isolated 
horizons [8]? In the spherically symmetric case, these issues are straightforward to address 
and an essentially complete solution is known. It would be very interesting to answer these 
questions in the axi-symmetric case. 



dM 



8ttG 



da + Vt Q dJ . 
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ii) In the analysis of section V, let us drop the restriction to t a Q and consider general 
permissible vector fields t a . Unlike the vector fields £?n = N r £ a , the vector field t a is not 
necessarily causal. Therefore the matter flux j AH T a i,t a T b d 3 V need not be positive. Similarly, 
if ip a is not a Killing field of qw,, the gravitational flux need not be positive. Therefore, 
although the area a always increases with R, E t can decrease as R increases. This is the 
analog of the Penrose process in which 'rotational energy' is extracted from the dynamical 
horizon. Do Einstein's equations with physically reasonable matter allow one to extract all 
the rotational energy? Once the first question in i) above is answered, one would have an 
essentially complete control on such issues in fully dynamical processes. 

iii) In a gravitational collapse or a black hole merger, one expects the dynamical horizon 
in the distant future to asymptotically approach an isolated horizon. Is this expectation 
correct? If so, what can one say about the rate of approach? There exists a useful charac- 
terization of the Kerr isolated horizon [14]. Under what conditions is one guaranteed that 
the asymptotic isolated horizon is Kerr? On an isolated horizon one can define multipoles 
invariantly [15] and the definition can be carried over to each cross-section of the dynamical 
horizon. Can one physically justify this generalization? If so, what can one say about the 
rate of change of these multipoles? Can one, for example, gain insight into the maximum 
amount of energy that can be emitted in gravitational radiation, from the knowledge of the 
horizon quadrupole and its relation to the Kerr quadrupole? Is the quasi-normal ringing 
of the final black hole coded in the rate of change of the multipoles, as was suggested by 
heuristic considerations using early numerical simulations? 

iv) As the vast mathematical literature on black holes shows, the infinitesimal version 
(5.4) of the first law is conceptually very interesting. The finite balance equation (5.3) is 
likely to be even more directly useful in the analysis of astrophysical situations. In particu- 
lar, there exist an infinite number of balance equations. Can they provide useful checks on 
numerical simulations in the strong field regime? Similarly, the Hamiltonian framework of 
Booth and Fairhurst could be used as a point of departure for quantum mechanical treat- 
ments beyond equilibrium situations. Can one extend the non-perturbative quantization of 
[15, 16] to incorporate these dynamical situations? To naturally incorporate back reaction 
in the Hawking process? 

v) There has been considerable interest in the geometric analysis community in the in- 
equalities conjectured by Penrose, which say that the total (ADM) mass of space-time must 
be greater than the area of the apparent horizon on any Cauchy slice. In the time symmetric 
case (i.e., when the extrinsic curvature on the Cauchy slice vanishes) this conjecture was 
proved last year by Huisken and Ilmamen [18] and Bray [17]. The area law of dynamical 
horizons provides a nice setting to extend this analysis not only beyond the time symmetric 
case but to establish a stronger version relating the area of apparent horizons to the future 
limit of the Bondi energy at null infinity. 
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